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ABSTRACT 

The observed rapid onset of star formation in molecular clouds requires rapid formation of dense 
fragments which can collapse individually before being overtaken by global gravitationally-driven 
flows. Many previous investigations have suggested that supersonic turbulence produces the nec- 
essary fragmentation, without addressing however the source of this turbulence. Motivated by our 
previous (numerical) work on the flow-driven formation of molecular clouds, we investigate the ex- 
pected timescales of the dynamical and thermal instabilities leading to the rapid fragmentation of 
gas swept up by large-scale flows, and compare them with global gravitational collapse timescales. 
We identify parameter regimes in gas density, temperature and spatial scale within which a given 
instability will dominate. We find that dynamical instabilities disrupt large-scale coherent flows via 
generation of turbulence, while strong thermal fragmentation amplifies the resulting low-amplitude 
density perturbations, thus leading to small-scale, high-density fragments as seeds for local gravity to 
act upon. Global gravity dominates only on the largest scales; large-scale gravitationally-driven flows 
promote the formation of groups and clusters of stars formed by turbulence, thermal fragmentation, 
and rapid cooling. 

Subject headings: gravitation — instabilities — turbulence — methods:analytical — stars:formation 
-ISMxlouds 



1. MOTIVATION 

There is increasing evidence that star forma- 
tion (in the solar neighborho od) follows rapi d ly on 
formation (IHartmann et all 120011: 

3~ 



molecular cloud 

I2007t lElmegreenl 
2007). This evidence suggests that the density en- 
hancements in which stars form are produced during 
the cloud formation phase - and not after the cloud 
has formed, as implicitly assumed by the initial condi- 
tions of many numerical models. Global gravitational 
modes can sweep up mate rial at the edges of clouds 
(|Burkert fc Hartmannll2004f). leading to the formation of 
filam ents (jVazquez-Semadeni et alJ l2007t iHeitsch et al.1 
120081 ). Yet for local gravitational collapse to win, the 
clou d needs to be seeded with lo c al perturbations early 
on (Burkert fc Hartmannl l2004t IHartmann fc Burkert] 
2007). Not only need these density enhancements to 
arise early, but they also need to be non-linear, i.e. the 
cloud must acquire strong small-scale density perturba- 
tions during its formation. Thus, a rapid fragmentation 
mechanism other - and faster! - than gravity is needed. 

The scenario of flow-driven cloud formation, in which 
large-scale shocked flows of atomic hydrogen in the warm 
neutral medium (WNM) assemble molecular clouds 
(|Ballesteros-Paredes et "ail 119991; IHartmann et al.ll2001f ). 
offers an elegant mechanism for the rapid fragmen- 
tation of the proto-cloud and the generation of non- 
linear density seeds. The key to the rapid fragmen- 
tation is a combination of stron g dynamical and ther- 
mal i nsta bilities in shocked flows dHennebelle fc Peraultl 
1 19991 and IKovama fc Inutsukal l2000l fo r one-dimensional 
models; IKovama fc Inutsukal 120021 for 2D models 
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including H? formation : lAudit fc Hennebe Id 120051 . 
Heitsch et alj|2005l l200l iHennebelle fc Auditll2007l and 
Hennebelle et al.ll2007l emphasi zing the formation of cold 
neutral medium clouds in 2D; IVazauez^S cmadcni et ail 
20061 Vazquez - Semadeni et al.1 [2007 and Heitsc h et al.l 
20081 for three-dimensional models, the latter two includ- 
ing self-gravity). Rather than generating turbulence by 
imposing e.g. a velocity field chosen ad hoc in Fourier 
space, density and velocity structures arise naturally dur- 
ing the formation of the cloud in this scenario. 

In this study, we attempt to throw some light on the 
relative importance of various classes of instabilities at 
play, via a consideration of their timescales. We aim to 
assess whether fragmentation processes other than grav- 
ity are acting rapidly enough during the formation of 
molecular clouds to allow the rapid onset of local star 
formation. Therma l effects - s trong cooling and ther- 
mal instability (TI, iFieldl 119651 ) play a dominant role, 
weakening the effective equation of state in a shocked 
flow and thus allowing rapid fragmentation. The TI de- 
rives its strength from a combination with turbulence 
which generates density perturbations on small scales. 
The resulting rapid growth of small-scale perturbations 
allows the rapid onset of local gravitational collapse, be- 
fore global collapse modes overwhelm any pre-existing 
low-amplitude density variations. Hence, in the context 
of flow-driven cloud formation, thermal fragmentation is 
the key to the rapid onset of star formation. 

The expressions for the timescales and the instability 
conditions are given in £jH We discuss the parameter 
regimes for the instabilities and their consequences in £j3] 
and summarize in $4] 

2. TIMESCALES AND INSTABILITY CONDITIONS 

The instabilities enabling the rapid fragmentation 
of shocked flows are associated with characteristic 
timescales. These will serve as a vehicle to estimate the 
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relative importance of the instabilities for fragmentation 
and turbulence generation. 

Since we are mainly interested in the competition be- 
tween global gravity and local fragmentation processes, 
we restrict ourselves to the discussion of four instabili- 
ties which can be split into two groups of physical pro- 
cesses, namely condensation processes (driven by cool- 
ing [ 32.1] and gravity [ 32.2] ) and generic fluid insta- 
bilities (ram pressure imbalance [ 32.3| and shear flows 
[ 32.4] ). This list of instabilities is by no means exhaus- 
tive ( 32. 5j) . Our choice is to some exte nt guided by 
the r esults of our numerical simulations (jHeitsch et al.l 
2008), and by the notion of cloud formation in large- 
scale, organ ized flows such as spiral arms of ga laxies 
(iKim et al.ll2003t|Kim fc Ostrikeijl200l lElmegreenll2007l: 
iDobbs fc Bonnelll 120071. 12008f). or ex panding supernova 
sh ells (e.g. [El mcgrccn fc Ladd [19771 see also discussion 
in IBallesteros-Paredes fc Hartmannl 2007h . 



2.1. Thermal Instability (TI) and Cooling Timescales 

The TI rests on the balance between heating and cool- 
ing processes. The astrophysically most relevant mod e 
is the isobaric condensation mode (e.g. iBalbusI [l995f l. 
The condensation mode's linear growth rate is indepen- 
dent of the wave number as long as the perturbation 
can adjust to an isobaric state and under neglection 
of heat conduction, h owever, perturbations on smaller 
scales will grow first (|Burkert fc LinlfeOOOf) . Thus, in 
the linear regime, the condensation mode is limited to 
scales smaller than the sound cros sing length (see also 
Vazquez-Semadeni et al.l 120031 and iHennebelle fc Audit! 
20071) . 



where 
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is the cooling timescale with temperature T and particle 
density n. Thermal energy losses due to optically thin, 
collisionally excited atomic lines are given by the cooling 
function A(T) in erg cm 3 s . The sound speed is 
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with the mean molecular mass u. Since the TI spans 
a range of 2 orders of magnitude in temperature, con- 
dition (fT]) can become quite restrictive, but it varies 
strongly with temperature and the strength of the cool- 
ing (see Fig.[T|). The scale A c can drop to a few tenths of 
a parsec for parameters typical of the WNM. 

In the non-linear regime, and on scales A > A c , the gas 
can still cool isochorically such t hat the result ing pres- 
sure drop will generate waves (e.g. lBalbudll99l . leading 
to additional fragmentation. In that sense, equation (fT]) 
is not a strict upper limit for the TI, since it only refers 
to the (linear) evolution of the condensation mode. Nev- 
ertheless, we will use it as a proxy for estimating the 
importance of the TI since it gives the largest scale out 
of which a single coherent cold fragment can form and 
thus is emphasizing the local nature of the TI. On small 
scales , the condensation mode is limited by the Field 
(|1965f ) length, below which heat conduction will become 
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Fig. 1. — Composite cooling curve (see text) for various ioniza- 
tion degrees and for metallicities corresponding to the solar neigh- 
borhood, the LMC and the SMC. 



important (see e.g. Koyama & Inutsuka 2004). The con- 
densation mode will grow if 
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where £ nA — T is the loss-heat function in erg s 1 

(|Fieldlll965l) . 

The growth timescale of the condensation mode under 
the conditions |T]) and (j4]) is given by 
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-fk B T 
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(5) 



with the partial derivatives dC n = dC/dn and d£r = 
dC/dT. Equation ([5]) includes the condition for the con- 
densation mode, i.e. for tti < 0, the condensation mode 
does not grow. 

Thermal effects still may be important in the absence 
of the TI, in which case the thermal timescale is given 
by the more general cooling time (eq. [2]). A short r c 
does not necessarily entail TI. Rather, a short t c in the 
presence of a more or less constant heating term means 
that the gas cools rapidly to a minimum temperature and 
then stays at that temperature effectively isothermal, i.e. 
the effective adiabatic exponent 7 e // = 1, as in P oc 
rfleff. A j e ff < indicates thermal instability: with 
increasing density, the pressure drops. For < "f e ff < 1, 
fragmentation can still be enhanced in the presence of an 
external (gravitational or ram) pressure. 

We derive the cooli ng function from a combin ation 
of the rates quoted by iDalgarno fc McCravl (ll972T ) and 
IWolfire et all (fl995l) for T < 10 4 K while for T > 10 4 K 
we us e the tabulated curves of ISutherland fc Dopital 
([l993h . Figure Q] summarizes the cooling curves for 
a range of ionization degrees Xi a nd metallicities, cor- 
respo nding to solar abundances l|Dalgarno fc McCravl 
Il972f) and abu ndances typical for the LMC and SMC 
(|Garnetdll999h . 

2.2. Gravity, Jeans Instability ( GI) 

A (one-dimensional) region becomes gravitationally 
unstable if 

{c s kf 



P > 



AwG 
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with the wave number k and mass density p = \m. The 
corresponding time scale is 

T GI = {l-KGp-<? s k 2 )-V\ (7) 

again including the condition for instability. Gravita- 
tional collapse occurs for tqj € K. 

2.3. Non-Linear Thin Shell Instability (NTSI) 

The NTSI (|Vishniad 11991 arises in a shock-bounded 
slab, when ripples in a two-dimensional slab focus incom- 
ing shocked material and produce density fluctuations. 
The growth rate of the NTSI is mostly controlled by krj, 
the product of the wave number of the slab perturba- 
tion k — 2it/L, and the amplitude of the slab's initial 
displacement rj (equivalently, the amplitude of the colli- 
sion interface's geometrical perturbation). The instabil- 
ity is driven by lateral transport of longitudinal momen- 
tum, i.e. if the inflow is parallel to the x direction, and 
the slab is in the y-z-plane, i-momentum is transported 
laterally in y (and z), collecting at the focal points of 
the perturbed slab. The efficiency of lateral momentum 
transport is key to the development of the instability, 
since it is the imbalance of ram pressure at the focal 
points that eventually propels matter forward, d riving 
the growth of the slab's perturbation. IVishniacl (|1994f ) 
derived a growth timescale of 

tntsi w [c s k{kri) 1/2 ]- 1 - (8) 

iBlondin &: Marksl (|1996f) found that at constant rj and 
for small k, equation |8]) yields only a lower limit, while 
for large fc, the analytical growth rates agree well with 
the numerical results. The reason for this seems to lie in 
the efficiency of deflecting the incoming flow: for small k, 
a small fraction of the incoming flow's momentum is con- 
verted to lateral motions, while a large part compresses 
the slab (depending on the equation of state, this could 
lead to an increase in energy losses). For the parameter 
study (33j, we set rj = L/4. Keeping k constant and 
varying rj affects the NTSI in a similar way as keeping rj 
constant and changing k, with a weaker dependence on rj. 
Magnetic fields can suppress the NTSI if the ram pressure 
of the incoming flow drops below the magnetic pressure 
(|Heitsch et al.ll200l . The NTSI is an efficient dynamical 
focusing mechanism for large-scale gas streams, resulting 
in the rapid build-u p of massive cores in the focal points 
( Hueckstaedt 2003) as possible sites for massive star for- 
mation (jHeitsch et al.|[2008l) . 

2.4. Kelvin- Helmholtz- Instability (KHI) 

The NTSI converts the highly compressible modes of 
the inflows into shear flows at the flanks of the perturbed 
slab. These shear flows give rise to the KHI, which thus 
is secondary to the NTSI, but which also is the main 
turbulence generation mechanism. Th e turbulence i n 
turn leads to the saturation of the NTSI (|Vishniadll99ll ). 
In the simplest — inc ompressible - scenario of the KHI 
( Chandrasekhar 1961), the shear layer has constant den- 
sity, and the velocity profile is a step function. In this 
case, the growth timescale is given by the velocity differ- 
ence AU and the wave number as 

Tkhi = (kAU)~\ (9) 

i.e. the system is unconditionally unstable. This is defi- 
nitely the most extreme case: magnetic field components 



parallel to the shear flows can stabilize against the KHI, 
and for compressible flows, the system will be stable for 
all those wave numbers wh ose effective M ach number is 
larger than a critical value (lGerwinllT9 68). For a detailed 
discussion of the KHI, see lPalotti et all ()2008l ). 

2.5. Other Instabilities 

The list of instabilities and their effects con sidered here 
is by no means exhaustive. For instance, llnoue et alJ 
(2006) discuss a corrugation instability of evaporation 
fronts appearing at the interface between the WNM and 
the CNM on timescales of approximately ~ 1 Myr in the 
WNM, and down to 0.01 Myr in the CNM. This - so they 
argue - is on the order of the cooling timescales in each of 
the media, so that the instability could contribute to the 
generation of turbulence in the ISM. Hence the corruga- 
tion instability will play an important role for the dynam- 
ics of the WNM /CNM Howe ver, as equation (47) and 
Figure 9 of llnoue et al.1 (|2006f) demonstrate, the growth 
timescale of the corrugation instability is in the (for our 
application) interesting range of < lMyr only for sub- 
parsec scales, which are substantially smaller than the 
large-scale flows sweeping up whole clouds as envisaged 
in the current study. 

In a similar vein, we do not discuss turbulence 
generation by the TI alone. Thi s possibility ha s been 
investigated exten si vely (e .g. iBurkert fc Linl 120001: 
Kritsuk fe Normanl l2002bl ; IVazquez-Semadeni et al.1 
2003?) . ICoxl (|2005f) points out that because of the 
limited energy reservoir, thermal pressure variations in a 
bistable ISM are expected to have negligible dynamical 
effects, requiring additio nal sources or triggers for 



'lggcrs to) 
nl l2002ar i. 



turbulence generation (e.g. iKritsuk fc Norma 



2.6. Criteria for Instability Dominance 

The characteristic timescales (eqs. |5I7I8I9| ) can be 

used to derive criteria for the dominance of a given insta- 
bility in terms of physical quantities. We chose (particle) 
density and temperature, to facilitate a straight-forward 
comparison to ISM regimes. Attributes in square brack- 
ets in each of the following subsection titles stand for 
the line-styles used in Figure O Not all combinations 
are physically relevant. For instance, r c <C tqi generally 
for gravitationally unstable regimes ( §3.2|) . The ratio of 
specific heat capacities is set to 7 = 5/3, and the mean 
molecular weight = 2.36. Obviously, these choices are 
not valid over the full parameter range discussed, but 
they will considerably simplify the discussion. 

2.6.1. Gravity threshold [green dashed] 
A region will be gravitationally unstable if 



n > 



lkB k 2 T 



inGji 2 



(10) 



- 7 - 4x10 m(r' cm 

This is just equation ^ expressed in terms of tempera- 
ture. 

2.6.2. GI dominates NTSI and KHI [red dashed] 
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If condition (fTOjl is fulfilled, gravitation will dominate 
the NTSI if 



3. DISCUSSION 



(ii) 



where r\ is the displacement of the slab. Note the strong 
dependence on the wave number k, which hints at the 
power of dynamical fragmentation to prevent global grav- 
itational collapse. 

The condition for the KHI is similar to conditions AT " 
and (fTTj) . The GI dominates over the KHI for 
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(12) 



L 
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where M. is the Mach number of the shear flow. This 
is the most beneficial case for the KHI: for supersonic 
flows, the KHI is not unconditionally unstable anymore. 

Since conditions pip and (TT2")) each depend linearly on 
the temperature, they will behave similarly in the (n, T), 
plane. Without loss of generality, we can combine them 



2.6.3. TI dominates NTSI and KHI [red solid] 
The TI dominates the NTSI for 

1/2 
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This just mirrors the fact that with increasing density 
the cooling becomes stronger. A similar condition can be 
derived for the KHI by replacing (fcry) 1 / 2 by the velocity 
difference multiplied by the wave number. 

2.6.4. Sound crossing time limit for the TI [blue solid] 

As discussed in M2. H we use the sound crossing scale 
(eq. [T]) for the TPs condensation mode as a proxy to 
estimate the importance of the TI. It leads to a density 
threshold above which the condensation mode of the TI 
cannot be excited: 



n < 



1/2 



(k B T) 



3/2 



A(T) 



(14) 



i 2.2 x 10" 
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T_\ 3/2 / 10 26 A(T) 
[KjJ \[crg cm 3 s- 1 ] 



The smaller the scales, the less restrictive is the limit on 
the density for the condensation mode. The temperature 
dependence is given by the detailed form of the cooling 
curve. 



The conditions summarized in £12.61 allow us to iden- 
tify the dominant instabilities in the two-dimensional 
(n, T) parameter plane. Most of the conditions (TlT))) - (|14p 
depend on the spatial scale. Figure [2| summarizes the 
regimes in four diagrams, at characteristic (fixed) scales 
of 0.1, 1.0, 10 and 10 2 pc. Table Q] provides a key to the 
line styles and colors used in Figure [21 We will discuss 
the dominant instabilities ("JOrGSI) and an evolutionary 
sequence of a fluid parcel in the ISM f £|3.5p . 



3.1. Thermally Dominated Regime 

The most prominent feature in Figure [2] is the large ex- 
tent of the thermally dominated parameter regime (right- 
slanted hashing; blue and black). This regime is lim- 
ited by the growing importance of dynamical instabilities 
(red) towards low densities, and by the Jeans condition 
(eq. [TO)) towards high densities. Towards low densities 
(red), the thermal timescales get longer than the dynam- 
ical timescales, so that the gas will tend to behave adi- 
abatically. Within the blue hashed region, cooling dom- 
inates dynamical instabilities, leading to the TI within 
the temperature ranges indicated by the thick black hor- 
izontal bars. Outside the strictly thermally unstable re- 
gions, cooling still dominates and can lead to fragmen- 
tation when an external (ram or gravitational) pressure 
is applied. Thus, while the TI dominates the dynam- 
ics only within certain temperature ranges, thermal ef- 
fects (strong cooling) continue to dominate all through 
the black hashed region at higher densities. The up- 
per edge of the blue ribbon is determined by the sound 
crossing time condition (|14p . which is - strictly - only 
applicable to the TI. There are two thermally unstable 
regimes (according to eq. U): The lower one between 
10 2 < T < 10 4 K connects the warm and cold neutral 
medium (WNM and CNM), while the upper (of lesser 
strength) in the range of T « 10 5 connects the hot and 
warm ionized medium (HIM and WIM). 

The curved thick black line between 10 < T < 10 4 K 
denotes the thermal equilibrium curve, where heating 
terms balance cooling terms. At T w 10 4 K and 
n < 1 cm~ 3 , the ISM behaves quasi-isothermal. The 
signature of the thermal instability is a pressure loss 
with increasing density (compare to dashed lines of con- 
stant pressure). Moving to high densities and lower tem- 
peratures, the effective equation of state tends back to 
isothermality but stays sub-isothermal, i.e. j e ff < 1. 
The thermal equilibrium curve is describing an approx- 
imate evolutionary sequence of a fluid element from the 
WNM to the CNM (see 

Before we compare the relative strengths of the insta- 
bilities, Figure [3] offers a more detailed view of the TI 
in the WNM. As orientation and for comparison with 
Figure [21 the lines of constant pressure and the sound 
crossing time limit (eq. |14j ) for the TI have been indi- 
cated in the same style as in Figure [2] The color-shaded 
region denotes the growth timescale (eq. [5]) of the TI 
(see color bar to the right of the diagrams) in all loca- 
tions where the TI can be excited. There are two points 
to notice here: 

(1) Looking at the upper left panel (L = 0.1 pc), 
the timescale tti can become extremely short; however, 
around the thermally unstable region, it is on the or- 
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Fig. 2. — Summary of the instability regimes in (n, T)-space, for four characteristic spatial scales as indicated in the top left corners. 
The diagonal black long-dashed lines denote lines of constant pressure as indicated. The thick black line between 1 < log T < 4 represents 
the thermal equilibrium curve due to a balance between heating and cooling. The horizontal thick black lines at the bottom of each panel 
denote the ranges of the two thermally unstable regimes. Fragmentation by cooling dominates in all regions with right-slanted hashing (red, 
blue and black), while gravity dominates fragmentation in the left-slanted hashed region (green). Red stands for dominance of dynamical 
instabilities over cooling/TI. Blue/black indicates where fragmentation by cooling dominates. The blue ribbon and the black horizontal 
lines together indicate ranges of strict thermal instability. See Table [T] for a summary of the line styles and colors. 



TABLE 1 
Key to Figure [2] 



line style 


meaning 


equation 


hashed, slanted right 


fragmentation by cooling or TI 




hashed, slanted left 


fragmentation by gravity 




red hashed 


dynamics dominate cooling 




blue hashed 


cooling/TI dominates dynamics 




black hashed 


cooling (but not TI) dominates dynamics 




green dash 


Jeans criterion 




red dash 


GI dominates NTSI 




red solid 


TI dominates NTSI 


H 


blue solid 


TI crossing time 





der of 1 Myr (light green, traced by the thermal equilib- 
rium curve), comparable to dynamical timescales. This 
leads to a noticeable amoun t of thermally unstable gas 
(e.g. iHeiles fc Trola nd 2003 ) for observations of neutral 
hydrogen. and | Gazol et al.ll200ll ISanchez-Salcedo et al.l 
l20f)l lAudTt fc Hennebellell200.1 iHeitsch et a.1.1 120061 for 
evidence from numerical models). For the isothermal 
branch at T « 10 4 K, the TI is essentially absent. Thus, 
the lower the density on the (warm) isothermal branch, 



the higher the compression due to e.g. shocks needs to 
be to trigger the TI. Since the post-shock density in an 
isothermal gas scales with the square of the Mach num- 
ber, this condition is not overly restrictive. 

(2) With increasing spatial scale (solid blue line in the 
remaining three panels of Fig. [3]), the upper density limit 
for the TI drops due to the sound crossing time condi- 
tion (TTJJ. Strictly, the TI will only be excited if there 
is a color-shaded region below the blue line, so to speak. 
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Fig. 3. — A more detailed look at the TI in (ra, T)-space, for four characteristic spatial scales as indicated in the top left corners. The 
diagonal black long-dashed lines denote lines of constant pressure as indicated. The thick black line between 1 < logT < 4 represents the 
thermal equilibrium curve due to a balance between heating and cooling. The horizontal thick black lines at the bottom of each panel 
denote th e ra nges of the two thermally unstable regimes. The blue line denotes the maximum density given by the sound crossing time 
limit (eq. 1141 ). The color-shaded region shows the growth timescale of the TI (eq. [5]). Light green corresponds to a timescale of ?s 1 Myr. 



This leaves us with the somewhat puzzling result that 
for scales L > 10 pc, the TI seemingly cannot be excited. 
The solution is two-fold: (a) At large scales (and lower 
densities/higher temperatures), the upper TI at 10 5 K 
kicks in, and (b) the notion of fixed scales is misleading. 
The TI will be triggered by compressions and/or tur- 
bulent mixing equivalent to a reduction of spatial scale. 
Moreover, since it is a condensation mode, we cannot re- 
gard the evolution of the TI at a fixed scale, but need to 
follow a fluid parcel (see §3.5| . 

3.2. Gravitation against Thermal Instability 

The Jeans criterion (eq. [10]) is fulfilled for all pairs 
(n,T) above the green dashed line in Figure [2]- e.g., for 
densities n > 10 5 cm -3 and temperatures T w 10 K at 
0.1 pc. With increasing scale, the (n, T)-regime dom- 
inated by gravitation extends more and more to lower 
densities and higher temperatures. At large scales (« 
10 pc), gravity starts to dominate the parameter space, 
indicating the importance of global gravitational modes. 

Figure 0] compares the relative strengths of gravita- 
tional over thermal condensation 3 . Shown are the same 

3 We are using "thermal condensation" for condensation due to 



plots as in the previous figures, but now overlaid with 
the ratio t c /tgi for the gravitationally dominated (n, T)- 
regime (eq. [7|). Yellow (see color bar of t c /tgi) de- 
notes r c = tgi- The cooling timescales are substantially 
shorter than any gravitational timescales over the whole 
parameter range accessible to the GI. If the TI can be 
excited (e.g between 80 < T < 5000 K), this means that 
thermal condensations will be the dominant fragmenta- 
tion process. For parameters (n, T) not subjected to the 
TI, the short cooling timescales mean only that any ex- 
cess energy due to compressions (by gravity!) will be 
efficiently radiated away. 

3.3. Gravitational against dynamical effects 

The relative importance of gravity and dynamics can 
be read off equations (fT0)) - (fT2|) . The discerning line for 
GI against NTSI is shown in Figure Has a red dashed 
line - obviously, the NTSI ceases to be important before 
densities high enough for gravitational collapse can be 
reached. Equations (fTQ|) and (|T2|) on the other hand show 
that for Mach numbers Ai > 1, the KHI can still dom- 
inate in gravitationally unstable gas. This is the regime 

strong cooling, not — as it occasionally happens in the literature - 
for gravitational collapse. 




Fig. 4. — Comparison of characteristic times for self-gravity and cooling in (n,T)-space, for four characteristic spatial scales as indicated 
in the top left corners. The diagonal black long-dashed lines denote lines of constant pressure as indicated. The thick black line between 
1 < log T < 4 represents the thermal equilibrium curve due to a balance between heating and cooling. The horizontal thick black lines at 
the bottom of each panel denote t he r anges of the two thermally unstable regimes. The blue line denotes the maximum density given by 
the sound crossing time limit (eq. 1141 1. The color-shaded triangles show the ratio t c /tqj (eqs. [5J, [7]) for the pairs of (n,T) subject to 
the GI. Yellow would correspond to t c = tqi - obviously, t c < tqi for all gravitationally dominated pairs of (n, T). 



of turbulent fragment ation in a Jeans-uns t able m edium 
discussed in detail bv lMac Low fc Klessenl ([2004) . 

3.4. Thermal against dynamical effects 

Dynamical instabilities dominate below the red (NTSI) 
and the black (KHI) lines (Fig. [2]), generally however 
in the red hashed region. They impose a lower density 
threshold on the importance of cooling (and of the TI 
within the applicable temperature range), i.e. the gas 
will tend to behave more and more adiabatically for those 
regions in (n, T)-space. 

Figure [5] shows the ratio of the cooling time scale 
(eq. [2]) over tntsi ( e q- [13] ) , where the latter is used 
as a proxy of dynamical instabilities. Since the cooling 
timescale does not depend on the spatial scale, dynam- 
ical instabilities can become dominant for small scales 
[L = 0.1 •••1.0 pc) in the WNM (T w 10 4 K). With 
increasing density, cooling dominates dynamics, but the 
TI is still limited by condition (T3j). 

Dynamical instabilities entail the generation of turbu- 
lence - especially at the high Reynolds numbers in the 
ISM. The turbulent cascade leads to a spectrum of den- 
sity perturbations, populating a wide range of spatial 



scales. In other words: instead of proceeding at a fixed 
scale, the TI will act on all density perturbations within 
its parameter range (see Fig. [3]) at scales smaller than 
the sound crossing scale simultaneously. 

3.5. An Evolutionary Sequence 

The strength and importance of the TI is based on the 
one hand on the fact that it grows on the smallest scales 
first (jBurkert fc Linll2000D and thus can feed on ubiqui- 
tous small-scale structures generated by turbulence. On 
the other hand, the condensation mode of the TI does 
not proceed at a fixed scale: a fluid parcel which gets 
thermally unstable will contract to smaller scales and 
thus "delay" running into the sound crossing scale con- 
dition (fT4"|) . 

To see this, consider the thermal equilibrium curve 
(solid black curve in e.g. Figs. [2] and [3]). Numerical sim- 
ulations of the fragmentation of shocked flows in a ther- 
mally unstable medium demonstrate that the gas gener- 
ally follows this thermal e quilibrium curve (see Fig. 3 of 
Audit fc Hennebelld 120051 and Fig. 14 of iHeitsch et all 
2006f ). so that it can be interpreted as an evolutionary 
sequence from the WNM to the CNM. While the equi- 
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Fig. 5. — Comparison of characteristic times for dynamical instabilities and cooling in (n,T)-space, for four characteristic spatial scales 
as indicated in the top left corners. The diagonal black long-dashed lines denote lines of constant pressure as indicated. The thick black 
line between 1 < log T < 4 represents the thermal equilibrium curve due to a balance between heating and cooling. The horizontal thick 
black lines at the bottom of each panel denote the ran ges of the two thermally unstable regimes. The blue line de not es the maximum 
density given by the sound crossing time limit (eq. 1141 ). The color-shaded regions show the ratio t c /tjjj (eqs. [2], 1131 ). Note that the 
ratio decreases monotonically with increasing density. 

also reduced the scale, so that effectively the thermally 
unstable regime (blue hashed) shifts upwards to higher 
densities, moving us from the 1 pc panel to the 0.1 pc 
panel in Figures [5] or [31 

Once the fluid parcel has entered the thermally unsta- 
ble regime, it evolves (more or less) along lines of con- 
stant pressure towards higher densities and lower tem- 
peratures. In the condensation mode, this is equivalent 
to a further shrinking of the fluid parcel's volume, so that 
the blue ribbon of thermal instability essentially "moves" 
to smaller scales with the fluid parcel on its way along 
the equilibrium curve. This "co-moving" or Lagrangian 
behavior is crucial for the strength and the importance 
of the condensation mode of the TI in the ISM. 

Yet the above statement calls for a caveat. A closer 
inspection of equation (fT4| reveals the following. As- 
suming that the fluid parcel evolves under constant 
pressure Pq - actually, for dynamical times substan- 
tially longer than thermal timescales, the fluid parcel 
will follow approximately the thermal equilibrium curve 
( Sanc hez-Salcedo et al.|[2002t and Fig. [J), which is sub- 
isobaric as indicated by 7e// < in the thermally un- 
stable regime - and with a constant cooling strength Ao, 



librium curve is well defined at the high-density end due 
to extremely short cooling time scales (see Fig. [3J , the 
scatter especially in the thermally unstable regime can 
be substantial and depends largely on the ratio between 
the d ynamical and the c ooling timescales (Figs. 14 and 
15 of lrleitsch et~alll2006[ ). 

A fluid parcel 4 in the WNM/WIM starts out at a scale 
of e.g. ~ 1 pc and a density of a little under 1 cm~ 3 . The 
parcel may be compressed due to a shock wave, or tur- 
bulence may lead to density variations. In any case, the 
density increases, propelling the parcel upwards along 
the isothermal branch of the equilibrium curve, reduc- 
ing the cooling timescales (Fig. [3J and thus eventually 
lowering the temperature. The fluid parcel will follow 
(approximately) the equilibrium curve along its sharp 
left-turn into the thermally unstable regime. Looking 
at the 1 pc panel, by then, the fluid parcel will have left 
the region of thermal instability (Fig. [2]). However, the 
compression above not only increased the density, but 

4 Strictly, the notion of a fluid parcel is slightly inappropriate. 
As the numerical models consistently show, the density structures 
tend to be filamentary rather than round. In that sense, the scale 
refers to the shortest axis of the filament/sheet. 
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the sound crossing scale condition ([14]) translates into 
an expression for the outer length scale Lt of thermal 
instability, 

L T <a n- 5/2 , (15) 
equivalent to equation (:1|, with the constant 



do = 2ir 



, \ 1/2 D 3/2 

lk B \ 1 P n ' 



M J A 

A spherical blob of constant mass mo will contract as 



(16) 



L v = b n- 1 / 3 , 



with the constant 



Thus, 



b 



6 mo 

7T /I 



1/3 



-13/6 



(17) 



(18) 



(19) 



meaning that the sound crossing scale (eq. pQ) shrinks 
more rapidly by a factor of w n 2 than the size of the fluid 
parcel. Even a one-dimensional compression Ly oc nT x 
would still result in Lt oc LynT^I" 1 . This would move 
the fluid parcel out of the thermally unstable regime. If 
the fluid parcel contracts "faster" than the sound cross- 
ing time criterion allows, the condensation mode will be 
suppressed, but the isochoric mode of the TI still will be 
excited. This is specifically the case for the black hashed 
regions (Fig. ^ where thermal timescales are substan- 
tially shorter than the dynamical timescales. The ther- 
mal pressure of the fluid parcel can then drop below that 
of the surrounding medium, and the resulting waves will 
lead to further fragmentation (smaller scales) and com- 
pression, triggering again the condensation mode. In 
that sense, while the co-moving picture is a simplifica- 
tion, the more realistic evolution of a fluid parcel involv- 
ing isobaric and isochoric modes will lead to even more 
fragmentation. 

Eventually, our fluid parcel will reach the sub- 
isothermal branch of the thermal equilibrium curve below 
« 80 K, entering the gravitationally dominated regime 
(as discussed earlier, this regime is characterized also by 
extremely short cooling timescales). In the presence of 
an "external" pressure (due to e.g. gravity), the sub- 
isothermal equation of state can lead to further fragmen- 
tation. 

However, gravity will not only play a local role. 
At large scales, global gravity can be dominant 
(Fig. [2 e.g. 10 pc). Global gravitational modes 
in combination with a finite extent of the cloud 
can easily lead to t he sweep-up of material a t the 
edges of the cloud (iBurkert fc Hartmanr] 120041: see 



IVazquez-Semadeni et al . 2007 and H eitsch et al, . 200lfor 
numerical evidence). Due to the rapid evolution of the 
local (thermal and dynamical) instabilities, the resulting 
filaments will have substructure, leading to further frag- 
mentation into cores. 

The thermal instability thus offers a short-cut to reach 
the gravitationally dominated regime from the warm, 
isothermal branch (at T = 10 4 K). Instead of having 
to compress the gas isothermally by a factor of ~ 10 4 - 
corresponding to a Mach number of 10 2 - only a modest 
compression is needed in the warm isothermal branch to 



move the fluid parcel up the thermally stable equilibrium 
line and "push it over the ledge" . This allows star forma- 
tion in environments without high pressurization and/or 
deep gravitational potentials. 

To summarize, the generation of turbulence due to dy- 
namical instabilities (red hashed regime) introduces den- 
sity variations at smaller scales, triggering the thermal 
instability, which in turn amplifies the density perturba- 
tions. The growth of the thermal instability is eventu- 
ally limited by the m ass reservoir available, and thus by 
the external pressure (|Kritsuk fc Normanll2002bl lal). This 
could take the form of global gravity, cloud collisions, ex- 
panding shells or galaxy mergers. Gravity acts locally 
on the high-density seeds, allowin g rapid collapse be- 
fore e dge effects of global gravity (|Burkert fc Hartmannl 
12001 overrun the local perturbations. 

3.6. Limitations 

Obvious limitations of the study include neglecting the 
effect of structure in the inflows, a generic cooling curve, 
the missing time-dependence and neglecting magnetic 
fields. Each of these is discussed in turn below. 

3.6.1. The Structure-less Inflows 

Our analysis approaches the problem of structured 
cloud formation from the angle of a "worst case" sce- 
nario, using the most unfavorable initial conditions to 
generate turbulence and substructure, by considering (as 
in our numerical work) the large-scale flo ws to be uni- 
form Obviously - as pointed out earlier (jHeitsch et alJ 
2006) - this assumption neglects the possibility that 
the flows themselves already might contain substruc- 
ture. This however is not problematic, as our goal is 
to demonstrate that substructure can arise from uni- 
form flows with only small, long-wavelength initial per- 
turbations. The introduction of structure in the in- 
flows - specifically a mixture of WNM/CNM instead of 
just WNM - would make it even easier to form sub- 
structure. Moreover, introducing substructure in the in- 
flows raises the question of just what form should be as- 
sumed. The most common initial conditions for numeri- 
cal simulations o f cloud formation com prise density per- 
turb ations fe.g. [inutsuka et all I2005D . velocity "noise" 
(e.g. IVazquez-Semaden et a l. 2007) , or a turb ulent ve- 
locity distribution (e.g. iHennebelle et al.ll2007l ). but the 
theoretical justification of these assumptions is unclear. 
Moreover, it is far from obvious that driving mechanisms 
such as that of a spherical stellar wind, H II region ex- 
pansion, or supernova bubble expansion are highly tur- 
bulent and structured in the absence of interaction with 
the surrounding medium. By ignoring substructure in 
the inflows we can consider the most general (and unfa- 
vorable) cases for generating turbulent substructure. 

3.6.2. Cooling Curve 

Our cooling curve at low temperatures is mainly lim- 
ited by the fact that it does not include molecular line 
cooling or the formation of molecules. Molecular cooling 
would lower the cooling timescales further, driving the 
cooling curve more to isothermality. The formation of H2 
lowers the sound speed by a factor of \[2 and the pres- 
sure by 2, thus introducing an additional pressure loss 
at high densities. While the thermal equilibrium curve 
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will describe the evolution of a fluid parcel approximately 
correctly from the WNM to the CNM, the details of the 
evolution at high densities and low temperature are less 
realistic. 

The heating and cooling functions T and A were de- 
fined for solar abundances and a Galactic UV back- 
ground. Reducing the metallicity by e.g. a factor of 10 to 
values typical for the SMC reduces the cooling strength 
below T « 10 4 K by approximately the same factor (see 
Fig- HI)- Thus, the blue ribbon in Figure would move 
up to higher densities. 

3.6.3. Time dependence 

Obviously, the analysis does not include any time evo- 
lution (although a rough idea of the dynamical evolu- 
tion can be gleaned from following the thermal equi- 
librium curve connecting the WNM and CNM, §3.5j) . 
Dynamical effects (turbulence, shocks) are expected 
to render thermal effects (cooling, TI) even stronger, 
specifically by providing small-scale seeds, and by al- 
lowing the TI to act on many scales simultaneously. 
This can clearly be seen in the high-reso l ution two- 
dimensional simulati o ns by Hcimcbc lle et alj (|2007t ) and 
iHennebelle fc Audit] (|2007f ). For a detailed discussion 
of (fully develop ed) turbulence in combi n ation with 
the TI, see e.g. IVazquez-Semadeni et al.l (|2003l) and 
I Audit fc Hennebelld (|2005l ). The latter authors point 
out that increasing the level of turbulence will lower the 
CNM fraction, an effect depending on th e ratio of dy- 
nami cal over thermal timescales (see also iHeitsch et al.1 
120061) . 

Certain simplifications had to be made, e.g. the Mach 
number for the KHI (eq. [9]) was set to 1. While this is 
a good approximation for many regimes in the ISM that 
exhibit trans-sonic flows, it will fail in cases of extreme 
shear flow conditions s uch as HVCs or mo lecular clouds 
(however, see Fig. 9 of IHeitsch et~aTll2006l) . 

3.6.4. Magnetic Fields 

Magnetic fields are not included. In relevance to the 
flow-driven cloud formation scenario, th ey would intro- 
duce a threshold crit erion for the KHI ( Chandra sekharl 
119611 ) and the NTSI (|Heitsch et alJl20q7ft. an d they may 
affect the evolution of the TI (TFieldJ |1965|). Critical 
field strengths and angles for the suppression of the 
condensation mode in one -dimensional converging flow s 
have been identified by (jHennebelle fc Perau lt 2000). 
Two-dimensional models of flow-driven cloud forma- 
tion with the magnetic field perpendicular to the in- 
flows (|Inoue fe Inutsukal l2008| ) demonstrate that while 
the fields can suppress the formation of high-density 
clouds, thermal fragmentation still occurs, leading to 
highly filamentary cold HI clouds. These results lighten 
somewhat the restrictive low limits on the perpendic- 
ular fie ld strengths derive d from one-dimensional argu- 
ments (jBergin et al.l 120041 ). On a larger (ga l actic s piral 
arm) scale, numerical models bv iKim et al.l (|2003| ) and 
iKim fc Ostrikerl (|2006h demonstrate the importance of 
magnetic fields for channeling gas streams. 

4. SUMMARY 

To form stars, molecular clouds need to fragment on 
scales vastly smaller than their overall dimensions. Self- 
gravity in combination with supersonic turbulence al- 



lows the rapid fragmentati o n of J e ans-unstable regions 
(iLarsonl 119811 : lElmegreenl 119931 : iMac Low fc Klessenl 
I2004T ) . Yet unless the density perturbations in molec- 
ular clouds are non-linear before the onset of global 
gravitational collapse, global gra vity tends to win and 
"over run" any local perturbation (jBurkert fc Hartmannl 
l2004f h Thus, a strong fragmentation mechanism is 
needed during the formation of the cloud to provide si- 
multaneously the observed turbulence as well as the ini- 
tial density seeds for the rapid local gravitational collapse 
mandated by the observed rapid onset o f star formation 
(jHartmann et al.ll200lT : lHartmannll2002f ). 

The flow-driven formation of molecular clouds in large- 
scale flows of atomic hydrogen provides a natural frag- 
mentation mechanism due to a combination of dynamical 
and strong thermal instabilities (see references in JQ) . As 
demonstrated by numerical models, this fragmentation 
mechanism does not need to recur on imposed turbulent 
velocity fields, but arises from the physical conditions in 
the WNM/CNM. The rapid fragmentation of the shocked 
flows will form high-density seeds for local gravity to take 
over, while gravitational forces o n the cloud scale c an 
sweep up material into filaments (IHeitsch et al.ll200ll ). 

The above outline of the connection between cloud 
formation and the initial conditions for star formation 
mainly rests on evidence from numerical simulations. 
In the present study we discussed the timescales of the 
instabilities dominating the evolution of shocked flows. 
Guided by our earlier work, we identify four instabili- 
ties, namely two condensation mechanisms (gravity and 
thermal instability) and two fluid instabilities (NTSI and 
KHI). We determine the parameter regime in density, 
temperature and scale in which each of the instabilities 
is expected to dominate. These are our findings and their 
implications: 

1. The (dynamical and thermal) instabilities leading 
to the rapid fragmentation of shocked flows dom- 
inate on small scales for reasonable parameters of 
the WNM and CNM (Fig. HJ). 

2. Cooling timescales (of ps 1 Myr in the thermally 
unstable regime) are generally substantially shorter 
than gravitational timescales for that same density 
and temperature regime (Fig. [4]). Thus, thermal 
fragmentation dominates gravitational fragmenta- 
tion during cloud formation. 

3. The strength of the condensation mode of the TI 
derives from two effects, namely that it grows on 
the smallest scales first, and that it is essentially 
a co-moving (or Lagrangian) mode, i.e. the un- 
stable region shrinks, thus keeping "longer" below 
the sound crossing scale A c (eq. [H]). Feeding on 
turbulent density perturbations below A c will allow 
the TI to grow on many scales simultaneously. 

4. The preference of small scales by the TI entails the 
preference to form sheets and filaments: the shorter 
axes are more susceptible to become unstable. This 
is an obvious parallel to the GI, with the - equally 
obvious but crucial - difference that the GI will act 
globally and locally. 

5. The fragmentation due to thermal effects is crucial 
for the formation of small-scale, high-density per- 
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turbations during the formation of the molecular 
cloud, to provide the seeds for rapid local gravita- 
tional collapse before global gravity dominates the 
dynamics completely. 

6. The all-important role of the TI due to the in- 
termittent nature of the turbulence driven by dy- 
namical instabilities allows an early fragmenta- 
tion of the converging flows, before (local) grav- 
ity can take over. This supports numerical find- 
ings that the (molecular) core mass spectr um might 
be set to some extent ()Dib et al.l [2008) early on 



during cloud formation by thermal fragme ntation 
(jHennebelle et al.ll2007HHeitsch et al.ll2008l ). 
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